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We consider the environment-affected dynamics of N self-interacting par- 
ticles living in one-dimensional double wells. Two topics are dealt with. First, 
we consider the production of entangled states of two-level systems. We 
show that by adiabatically varying the well biases we may dynamically gen- 
erate maximally entangled states, starting from initially unentangled product 
states. Entanglement degradation due to a common type of environmen- 
tal influence is then computed by solving a master equation. However, we 
also demonstrate that entanglement production is unaffected if the system- 
environment coupling is of the type that induces "motional narrowing". As 
our second but related topic, we construct a different master equation that 
seamlessly merges error protection/detection dynamics for quantum informa- 
tion with the environmental couplings responsible for producing the errors in 
the first place. Adiabatic avoided crossing schemes are used in both topics. 
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I. INTRODUCTION 

Real quantum systems are always coupled to environments. Research in the last two or 
three decades has shown that system-environment interactions greatly affect quantal coher- 
ences in the system. Open-system dynamics can often be described by a master equation 
for the reduced density matrix of the system. The full Hamiltonian is written as the sum 
of system-only, environment-only and system-environment interaction terms. The master 
equation for the system then depends on the system-only Hamiltonian plus other terms de- 
rived from the fundamental system-environment couplings. These additional pieces induce 
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an effective non-unitary evolution for the system that can, for example, lead to an apparent 
loss of quantal coherence. 

In many studies, the system-only Hamiltonian is taken to be the sum of one-body terms. 
An example of relevance to this paper arises for a system of N particles, each in a one- 
dimensional symmetric double well. The system Hamiltonian 



describes particle i oscillating with frequency U{ between the left and right sides of its well 
independently of the other particles in the system. When system-environment coupling is 
added, this behaviour can be drastically modified, with the quantum Zeno |jj or freezing 
phenomenon being the most extreme example. 

The purpose of this paper is to contribute to a study of closed and open quantum system 
evolution that goes beyond one-body system Hamiltonians by including interaction terms 
between the system particles. This is a vast and extremely rich field of enquiry, with many 
possible lines of development. For reasons to be explained below we will focus on two topics. 
The first is entanglement creation through adiabatic evolution and its degradation due to 
environmental influences. The second is quantum information error protection through 
adiabatic evolution. 

Common to both topics is the use of adiabatic avoided level crossing dynamics. This spe- 
cial case of quantum mechanical time evolution has applications in diverse areas of physics, 
from the Mikheyev-Smirnov-Wolfenstein effect for neutrino oscillations @ to quantum com- 
putation p. We use it here because it allows non-trivial dynamics to occur while retaining 
both an element of theoretical simplicity plus visualisability with the aid of level crossing 
diagrams. 

Whenever one has more than a single particle in one's quantum system, entanglement 
adds spice to the analysis. Entangled states are a unique feature of quantum mechanics, first 
studied in connection with the issue of non- locality ||. They are also of critical importance in 
discussions of decoherence, measurement and the quantum-to-classical transition ||. More 
recently, entanglement has emerged as a useful resource in the field of quantum information 
theory, enabling processes such as quantum cryptography || and teleportation |7j]. The 
production and manipulation of entangled states is a key element in any realisation of a 
quantum computer [||. 

As our first topic, we describe an efficient method for the preparation of entangled states 
which utilises the level crossing structure of the system's energy eigenstates. Specifically, we 
discuss the Bell states in the case of a bipartite system, W states in a tripartite system and 
generalised Wn states in multi-partite systems. We then subject the evolving systems to 
environmental interactions that affect the entanglement production process. We show that 
the outcome depends on the form of system-environment coupling: entanglement degrada- 
tion is the generic effect, but the special case of "motional narrowing" noise actually leads 
to zero degradation. Calculations presented in Ref. jnj simply took initial Bell states and 
computed the evolution of entanglement in the presence of environmental noise. The present 
work is a natural extension of this. 

As a separate but related topic, we show that a method of error protection using entangled 
states emerges from a variation of the above procedure. We devise a scheme to protect 



N 




(1) 



2 



the information encoded in a single qubit from bitflip and phaseflip errors. Our goal is 
to construct a master equation which seamlessly joins the dynamics of qubit encoding, 
error protection and decoding to the environment-coupling influences responsible for error 
generation. While the outcome of the open system evolution can be dissected into these 
quantum information theoretic subprocesses, the actual dynamics is simply what you get 
when you solve a certain master equation. 

The building block for our multi-partite systems is the two-state system of a particle in 
a one-dimensional double well, chosen for its generic properties and broad interest (see Ref. 
||). The well is taken to be partitioned by a barrier of height V that is large compared 
to the ground state energy E . (The two lowest energy levels will be collectively called the 
"ground state". For a symmetric well with V ^> Eq they are nearly degenerate.) Adopting 
quantum information theoretic parlance, we denote by "qubits" the states occupying the 
two-dimensional Hilbert space spanned by the left and right sides of the well (in the ground 
state). The left and right eigenstates, denoted |1) and |0), correspond to the +1 and — 1 
eigenvalues of a z , respectively. For an isolated well, the Hamiltonian has the form 

E x = E 1 + ua x + fa z , (2) 

where the energy splitting u is essentially a tunneling rate between the two sides of the well, 
and the parameter / is the well bias. 



II. GENERATION OF ENTANGLED STATES 

A. Two qubit system 

We begin with a bipartite two-state system, the simplest system in which one may discuss 
entanglement. Our aim here is to produce one of the Bell states, 

I^H^Oio^ioi)), (3) 

dynamically from certain initially factorisable states. 

Let us denote the two-particle (reduced) density matrix of qubits 1 and 2 as 

P = Pl®p2, (4) 

and consider a Hamiltonian of the form presented in Eq. (|2|) for each of the two qubits, 
together with a third term which couples them: 

H 2 = [ua xl + f{t)a zl ] <g> 1 + 1 <g> [uja x2 + f(t)a z2 \ + Xcr zl <g> a z2 . (5) 

Observe that the Hamiltonian is invariant under the interchange of the two qubits. For 
simplicity, only the well biases / are allowed to vary with time, and we demand that uo < 0.1A. 
This last condition is necessary in order for the desired level crossing structure to arise. (It 
is implicitly assumed that only the lowest states in each well are populated. We will need 
to relax this condition later on when we couple the system to a hot bath.) 
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The four eigenstates separate into an antisymmetric singlet state \^~) for all val- 
ues of u), f and A, and a symmetric triplet composed of linear superpositions of the set 
S = {|11), 1 00) , which asymptotes to S for large |/|'s. Note that the permutation 

symmetry of the Hamiltonian ensures that the evolution of the singlet state is disconnected 
from that of the triplet states at all times so we are dealing with an effectively three-state 
system. In Fig. [I], we plot the energy level diagram as a function of /. 
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FIG. 1. Level crossing diagram for the two-qubit system. The energy and the variable / are in 
units of A, and the dotted line denotes the decoupled state \^~)- 



The states |11) and 1 00) experience avoided level crossings with |\I/ + ) at / = —A and 
/ = A respectively, labelled points A and C in Fig. [I]. Hence, an initial state \<fi(t = 0)) = |11) 
created at / <C —A will continue to reside in the lowest energy eigenstate provided that / 
varies sufficiently slowly with time, and be adiabatically transformed, across point A, to the 
entangled state 

m = T)) = '* +> " vWTxj |n> + 7f(7^) |00> + 0(w2) ' (6) 

where T corresponds to any point in the period of time during which / lies in the range 
—A < / < A. Further adiabatic evolution beyond point C will turn \4>(t)) once more into 
a factorisable state 1 00) . Note in Eq. (H) that the entangled state \<p(t = T)) inevitably 
contains some order u/[f + A, / — A] "contamination" from the states |11) and 1 00) , which 
cannot be arbitrarily minimised without jeopardising the adiabaticity of the transitions at 
points A and C (see later). 

Figure § illustrates the evolution of the entanglement for | </>(£)), and we remind the reader 
again that the Bell state \^~) cannot be generated in this manner because of symmetry 
requirements. In this figure, the degree of entanglement is defined to be the entropy of 
entanglement E which is the von Neumann entropy of the reduced density matrix obtained 
by tracing out either of the two qubits [fllfl , 



E = -Trpi log 2 pi = -Trp 2 log 2 p 2 , (7) 

where Pi{p2) is the density matrix obtained after tracing out particle 2 (1). The related 
concept, entanglement of formation Ef, to be used below, is a generalisation of the entropy 
of entanglement for a mixed state of two qubits fi~2|| . 
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FIG. 2. Entanglement of \<fi(t)), where \4>(0)) = |11), as a function of time. We take the 
parameters u = 0.05A and f(t) = (—2 + i/2000)A, for arbitrary A. The time, t, is given in units of 
1/A. 



The entanglement will stay approximately constant at E ~ 1 when / is in the range 
—A < / < A. Note that if the two qubits had different oscillations frequences, mixing 
between the almost degenerate states is not prohibited by permutation symmetry (or 
the lack thereof). In this case, the degree of entanglement will vary significantly in the 
period —A < / < A, attaining the maximal value of E ~ 1 only at the point / = 0. 

Observe that if we had started instead in the state 1 00) which is approximately the 
highest energy eigenstate, we could in principle obtain the entangled state ^(|11) + 1 00) ) 
at the point / = 0. However, an inspection of the Hamiltonian 
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rewritten in the symmetrised basis {|11), |^ + ), 1 00) , I 1 ? - )}, shows that the splitting between 
the relevant energy eigenstates for the 1 00) # |11) transition at / = is of order uj 2 /X, 
as opposed to \f2uj for the |11) # transition considered earlier. It follows that the 

resonance width at point B is necessarily some cu/X times narrower than those at points A 
and C. Their respective adiabaticity parameters, defined as 
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where Ak is the splitting between the energy eigenstates, "res" denotes evaluated at reso- 
nance, and 7 > 1 indicates an adiabatic process, must also differ correspondingly by some 
factor of u 2 /\ 2 . Explicitly, 



characterise the three resonances. 

Hence, in order to manufacture the state ^=(|11) + 1 00} ) , a slowly varying / together with 
control over its accuracy to within ±u; 2 / A are necessary. Naively, the severity of the precision 
requirement seems easily alleviated by enlarging u. The benefit of a wider resonance width, 
however, is compensated for by a concomitant increase in the contamination from the \^f + ) 
state which contributes at order to/X. By comparison, the production of |^ + ) from |11) 
requires only that we stop the evolution somewhere between —A < / < A, and is therefore 
a much simpler task. 

1. System- Environment coupling: cold bath 

In the presence of environmentally-induced noise, we must solve a master equation for 
the reduced density matrix of the system.[] The form of the master equation depends on 
whether or not the environment is hot enough to induce transitions to excited states of the 
double wells. We first consider coupling to a cold environment, so that each particle remains 
in the ground state of its double well. 

The evolution of p is governed by 

j tP (t) = -i[H, P (t)\ - r rclax [c, [C,p\] - r rclax [C, [C,p\], (11) 

where ( and (' are Hermitian matrices obtained from the system-environment interaction 
Hamiltonian. For the sake of the example, we will choose £ = |( 1 + 0-3) ® 1, £' = 
1 (gi |( 1 + (73), which means that the relaxation terms arise from a system-bath interaction 
that couples to the left states |1) but not to the right states |0).Q Commencing with a 
pure state |11) and a small relaxation rate r re i ax , we show in Fig. ^| the evolution of the 
entanglement with the course of time. A "decay" in E in the period —A < / < A is evident. 
For a sufficiently large r re i ax , the density matrix simply becomes proportional to the unit 
matrix, and the entanglement goes to zero at all times. 

2. System- Environment coupling: hot bath 

We now consider coupling to an environment hot enough that the system-bath interac- 
tions will cause particles in the wells to populate states higher than the ground state. In 
a previous publication we considered in detail the effects of noise coming from a thermal 
bath coupled to the full excited state spectrum of a double well |L(J. Among the cases 



Certain conditions must be satisfied for the dynamics of an open quantum system to be well- 
described by a master equation. We consider only such cases in this paper. 

2 The other extreme has the environment coupling equally strongly to the left and right sides of the 
wells, so that Q and (' are proportional to the identity, and the double-commutator terms vanish. 
We will return to this point shortly. 
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FIG. 3. Entanglement of formation as a function of time in the presence of noise. We take 
r re iax = 0.5 x 10~ 4 A and display the time, t, in units of 1/A. The other parameters the same as in 
Fig. | 



considered was the effect of "side-blind" noise on a number of properties of the system. For 
the side-blind case, the bath couples equally to the left states |1) and the right states |0).| 
If we begin with a single particle on the left in the ground state of the symmetrical well, 
the thermal excitation of states that have very different R-L oscillation frequencies replaces 
with a multiperiodic tangle the previous simple oscillation of the probability of finding the 
particle on the left at a later time. The effect of a small system-bath interaction rate is to 
gradually damp these oscillations, leaving a time independent 50% occupancy on either side. 
As the ratio of the system-bath interaction rate to the oscillation rate is increased, however, 
the motion becomes more organized, and in the limiting case where this ratio becomes very 
large,0 we retrieve sinusoidal oscillation with a frequency that is given by the thermal av- 
erage of the frequencies of the individual levels JTO] . The analogous phenomenon is called 
"motional narrowing" in NMR studies. 

We will now look at entanglement-formation dynamics in a two-particle system with 
excited states. The L-R oscillation properties for a particular well-level are given by Eq. 



In the absence of excitation to higher energy states, side-blind noise does nothing of note. 

4 Note that by "hot bath", we refer to the situation where the temperature of the bath is high 
enough to excite states other than the ground state. The coupling, however, which determines 
the rate of these transitions, is a separate parameter. The "large T limit" arises when the time 
between successive system-bath interactions is much smaller than the period of the L-R precession. 
Observe that this condition can actually be fulfilled for any value of the coupling (and importantly, 
for a coupling sufficiently weak that one may treat it perturbatively in the derivation of the master 
equation) by making the free precession period sufficiently large. In the context of a double well 
oscillator this is achieved very easily, simply by adjusting the central barrier between the wells. 
Also note that this "large T limit" is identical to the limit where one obtains the familiar Quantum 
Zeno Effect, if the coupling to the environment is through a 03 operator rather than the identity. 
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(p|), but the oscillation parameters are different for the excited well states. We shall refer to 
the variable distinguishing these states as the "vertical" variable, as opposed to the (L, R) 
degree of freedom which we will term the "horizontal" variable. We will take the coupling 
among the vertical levels to again be "side-blind" , that is, not dependent on the horizontal 
variables, because that is the most interesting case. We can capture the essential features 
of the problem by considering only two vertical levels, with equal transition rates, T, each 
to the other (the latter as in the case of a high temperature for the thermal bath), in which 
case we make the replacements in Eq. (|TTD , p(t) — > p(Ei, t), i = 1, 2 on the lefthand side and 



rWC,[C,p]] - r(£,p(E j ,t)-2 P (E i ,t)) (12) 

3 

on the righthand side, as well as specifying, within H, different oscillation parameters for 
the two energy states. 

Using this equation we now look at the results for the two particle-entanglement for the 
cases of zero, moderate, and large system-bath interaction rates. The results are shown in 
Fig. 4. We see substantial degradation of entanglement-formation in the case of moderate 
interaction rate. But in the case of large interaction rate, we see a pattern of entanglement 
creation that is similar to that in the case with no noise. In this limit, the system is made 
to change vertical state so rapidly that it can only respond to an average Hamiltonian, with 



quantal coherence fully maintained |]10fl . For these reasons, the entanglement-formation 
dynamics succeeds in this limit. 
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FIG. 4. Entanglement of formation, Ef as a function of time, in the presence of "side-blind" 
noise. The solid, dotted and dashed curved are for zero, moderate amd large system-bath in- 
teraction rates, T = 0, 1 and 1000 respectively. For simplicity, we took a system with only two 
vertical states, and the initial condition was chosen to be |11), with both particles in their ground 
states. We have taken the parameters u(E ) = 0.05, v(Ei) = 0.1, X(Ei,Ej) = 20^Juj(Ei)uj(Ej) 
and / = — 3 + t/500, where u, /, A and T are in units of an arbitrary energy parameter E novm , and 
the time in units of l/E norm . Note that the position of the level crossing differs for the solid and 
dashed curves as A has been taken to be dependent on the value of the vertical variable. 



For completeness, we mention that if the bath does not couple with equal strength to the 
two horizontal states |0) and |1), the behaviour is qualitatively the same as that obtained 
from Eq. (|TI|) for the cold bath, where ( ^ 1 . 
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B. Multi-Qubit system 



Consider now the W states, 



\W) 001 = ^=(1001) + |010) + |100)), 

|iy) 110 = -^(|iio) + |ioi) + |on)). 



(13) 



These may be considered maximally entangled 3-qubit states, in the sense that each has the 



maximum amount of bipartite entanglement if any one of the three qubits is traced out |T3 



It turns out that the production of a three-qubit W state from an initially factorisable 
state such as | 111) via adiabatic transitions is completely analogous to the generation of the 
symmetric Bell state from |11) considered previously. We now elaborate on this point. 

The Hamiltonian for our three-qubit system is a direct generalisation of that in Eq. (^): 



H 3 = u [<j x1 + a x2 + a x3 ] + f(t) [a zl + a zX + a z3 ] + A [a zl a z2 + o z2 o zZ + a z3 cr zl ] 



(14) 



where, for simplicity, we are dealing explicitly with the permutation invariant case. Decom- 
position of the tensor product 2®2cg>2 = 4©2©2 reveals that the permutation symmetric 
set {|111), |W / )no, |W)ooi, 1 000) } can be treated in isolation, with its evolution governed by 
the 4x4 block in the Hamiltonian 



/3/ + 3A y/3u 
V3u /-A 
2uj 
\ 



\ 

2co 
-/ — A V3u 
V3u -3/ + 3A I 



(15) 



written in the symmetrised basis. The level crossing diagram for this set of states is plotted 
in Fig.|[ 
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Three resonances mark the evolution of the lowest energy eigenstate from / < — 2A to 
/ > 2A. To zeroeth order in u, the identity of the state is adiabatically transformed as per 
the evolutionary sequence 

|H1) ^ \W) 110 ^ |WOooi ^ 1000), (16) 

where the conditions above the arrows denote the points at which the transitions take place. 
Commencing with the product state 1 1 1 1 } at / < — 2A, the two W states can be easily 
obtained respectively in the domains — 2A < / < and < / < 2A. 

Note in passing that the state |GHZ) = |111) + 1 000) coincides approximately with 
an eigenstate at / = and can in principle be accessed from 1 000) or |111) via the highest 
energy eigenstate. However, as a direct analogue of the 4g(|ll) + 1 00) ) state considered 
in the previous section, the production of the GHZ state by our method suffers the same 
problem of non-adiabaticity and again requires precise control over the value of /. 

The generalised form of the W state for N qubits is defined as 

\W) N = -±=\N-1,1), (17) 

where \N — 1, 1) denotes the totally symmetric state with N — 1 qubits in the state |0) and 
one qubit in the state |1). 

The residual bipartite entanglement between any two pairs of qubits, upon tracing the 
other qubits out, has been shown to have the maximum value attainable in an N qubit 
system Since \W)n states are symmetric with respect to permutation of the N qubits, 
they can always be made from |11 • • • 1) or 1 00 • • -0) in the manner discussed above. 



III. ERROR DETECTION AND CORRECTION 

An important problem in the study of quantum information theory is the protection of 
information from various types of errors. Any practical realisation of a quantum computer 
will see the physical systems, the states of which comprise the qubits, experience effective vi- 
olation of unitary evolution and error generation, both due to coupling to their environment. 
A number of error detection and protection schemes have been devised that encode a given 
qubit through entanglement with other otherwise redundant qubits |J |16|] . Our goal in this 
section is to construct a master equation that seamlessly merges error generation, protec- 
tion through entanglement and detection within a single dynamical framework. While the 
resulting evolution can be usefully understood through reduction to information theoretic 
subprocesses, at the level of pure physics the single equation that one solves encapsulates 
both the quantum processes over which the experimenter has no control, and the knobs the 
experimentalist turns to change the biases. 

In what follows, averted level crossing schemes similar to those described in the previ- 
ous section are used to encode a qubit in a way that it is protected from certain outside 
disturbances, following the steps: 

1. We consider an initial pure state given by = a\ll) + b\01) = [a\l) + 6|0)]|1). Our 
goal is to preserve the information encoded in the first qubit. 
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2. We take a noise term that is derived from an external bath that flips qubits, indepen- 
dently, for one or both particles, and we follow the density matrix for the resulting 
mixed state. 

3. We transform the state in the time interval t — to t — t e by adiabatically changing 
the coefficients of certain terms in the Hamiltonian; we refer to these functions as the 
bias functions. Next we leave the system in this state for a length of time th ("hang 
time") long compared to t e . 

4. To reclaim the initial state we adiabatically restore the biases to their initial values 
and then measure the value of the second, control, coordinate. If the measurement 
gives back the original state |1), then for the first qubit we obtain the state that is the 
projection of our evolved density matrix on that state. The coding algorithm is such 
that this operation reclaims the density matrix of the initial pure state, with small 
degradations coming from the small non-adiabaticity due to doing the encoding and 
decoding over finite time intervals, or with degradations from letting the hang time be 
so large that multiple errors are likely to have accumulated. If, on the other hand, we 
measure the value |0) for the control state, then, beginning from the projection of our 
evolved density matrix on the state |0) a simple, and implementable unitary operation 
restores the initial density matrix. 

We begin this process by describing the system in terms of the four time dependent 
instantaneous eigenstates of the Hamiltonian that describes the system without noise, \£j)t, 
where j — 1,4. Defining , |£ 3 -) = \Cj)t=o, we take initial values, 

l£i)o= |i>® |i> 

|&)o = |0)®|l> 
l6)o= |1>® |0> 

|&)o = |0> ® |0> (18) 

The initial data is to be encoded in a combination of |£i)o, and |^2)o- We define a time 
dependent projection operator on the subspace used for encoding, 

p e (t) = ia>t<£iit+ik> t <6i* (19) 

We note that P = P e (0) = 1 ® |1)(1|, the projection operator on the value unity for the 
second qubit. We choose the Hamiltonian 

H = 2 f(t) <T x ®a x + a g ®<T g + 4 ! *[l- f{t)} a z ® 1 , (20) 

where f(t) is a function that begins at zero at t — 0, rises to unity at t — t e , remains unity 
until t = t e + th, and then decreases back to zero at t — 2t e + 1^. At t — 0, the eigenstates of 
H are the vectors shown in Eq.flTSp. We define U to be the time evolution operator for the 
interval t — to t — t e , under the action of this Hamiltonian. If the encoding time is chosen 
sufficiently large that the evolution is adiabatic, then for our encoding states we have, 

^i)o = |6k=e^ 1 (|ll) + |00))/v / 2 

f/|6)o=|6k = e*(|10)-|01))/v / 2 (21) 
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where the phases <pi and 02 are irrelevant to the remainder of the discussion. These remain 
eigenstates until t = t e + th- 

We will show that this dynamics protects against "bit-flip" errors defined by choosing 



the coefficients (, etc. in an equation like |ll] to be from the list of error operators a x ® 1 



and 1 ® a x , plus "phase errors" induced by 1 <8> a z and o z <g> 1 . (It is easy to modify the 
scheme to protect against all four bit-flip errors a x , y (g> 1 and 1 (g) cr x ,y, but at the cost of 
losing the phase error protection.) The key to this protection is that the any error operator 
(, picked from the above set has all matrix elements equal to zero in the subspace spanned 
by |£i)t e an d \^2)t e - Q We can express this property using the projection operator, P e (t), 

( = P e (t) C[l - P e (t)} + [1 - P e {t)](P e (t) + [1 - P e (*)]C[l " P*®] (22) 

which holds throughout the time interval t e < t < t e + th, noting that P e (t) is independent 
of time during this interval. 

In what follows, for clarity of exposition, we shall leave the noise turned off during the 
short encoding-decoding times; however in the numerical example presented later, we retain 
the noise throughout the process. We begin with a density matrix, 

P(0) = N6>o + 6|6>o] [a*<6|o + 6*<6lo] , (23) 

and let it evolve in a unitary fashion to t — t e , 

p(t e ) = U- 1 p(0) U = [a\^) te + b\£ 2 ) te ] [a*<6L + &*(6ld + n.a. (24) 

where "n.a." stands for small terms coming from non-adiabaticity, which can induce hopping 
into the ^3) and ^4) states. Between t = t e and t = t e + th the time evolution is governed 
by the equation, 

j t p{t) = -i[H(t e ),p(t)} - T[C , [C , p(t)}) = -i[H(t e ), p(t)} - r {p(t) - (p(t)(} , (25) 



where the noise term is of the same form as given in Eq. (|TT|) , the constant T sets the level 
of the noise, and ( is one of the error matrices enumerated above, with ( 2 = 1. We shall be 
concerned with the changes to the density matrix that are first order in T. However we note 
that if we omit the last term in Eq . (|25|) , the solution, p 1; would be just a damped form of the 
solution without noise, which is just the evolution under the Hamiltonian, H(t), constant 
in the interval t e < t < t e + t h , 

Pl (t e + t) = P e (t e ) exp[-iH(Qt] U- 1 p(0) U exp[iH(t e )t] P e (Q exp[-r*]. (26) 

The right and left projection operators are unnecessary here, since nothing has evicted us 
from the encoding subspace, but they are useful for the next step. For the effects on p(t e +th) 
of the ( dependent term in Eq . (|25|) , which we denote as p 2 , we can calculate the change in 



5 There will of course be nonvanishing matrix elements from this subspace to the subspace spanned 
by |£ 3 )i e and |£ 4 ) ie . 
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p(t e + th) to first order in T by direct time integration after inserting the (T = 0) form of 
p{t) from Eq. (|26D into the last term on the RHS of Eq. (|25|) . For our purposes all we need to 
know about this term comes from the projection operator structures in Eq.(p2|) and Eq.([26|) 
which directly give, 

p 2 (t e + t h ) = [1 - P e (t e + t h )\ p A [1 - P e (t e + t h )\ + OiT 2 ) (27) 

where pa is the result of the time integral discussed above. The total correction is given by 
Pi + P2- 

Next we decode, by adiabatically changing the parameters back to their values at t — 0. 
This is a unitary process implemented by the time evolution operator U for the interval 
t = t e + th to t = 2t e + th- We note that UP e (t e + th)U = Po, since the unitary evolution 
processes never cause a transition from one of the states \£j)t to another in the adiabatic 
limit, and the relative phase factors from the time evolution cancel in the construction of 
the projection operators. Thus we obtain, 

p(2t e + t h ) = P Q U exp(-iHth) U~ l p(0) U exp(iHt h ) U' 1 (1 - Tt h ) P 

+ (1 - P ) UpaU' 1 (1 - P ) + 0(r 2 ) (28) 

In the |£i)o, ^2)0 subspace the unitary matrix U exp(-iHth) U^ 1 is diagonal. Thus if we 
measure the second (control) variable and obtain unity (or no change from the initial value) 
there can have been no bitflip in the state of the first (information) state. Then, if we wish 
we can fine-tune the hang time, t h , so that U exp(— iHth) U' 1 in the |£i)o, ^2)0 subspace is 
a phase factor times the unit matrix. In this case we have recaptured the initial state of the 
information bit perfectly, using implementable transformations that are independent of the 
information itself. 

The above process can be characterized as being a limited error correction, in that we 
have corrected the errors that could have been created in the part of the final density matrix 
in which the value of the control bit is unity. In the case that the measurement of the control 
bit gives zero (this probability is itself of order Tth), then we can work out the explicit form 
for the second term on the righthand side of Eq. (p8|), and find a unitary operation that 
converts it to a multiple of the initial density matrix. This operation depends on which of 
the four error operators is used for the matrix (. Thus if we extended Eq. (^5|) to include 
several error terms at once, then, since we would not know which was operative in a given 
case, we could not implement the correction. Corrections of all of the error operators require 
encoding protocols involving more qubits, as discussed in Ref. [0, and in much subsequent 
literature. 

However, if we had replaced the double-commutator on the righthand side of Eq. (|25|) 
by a sum of terms, one for each of our protectable error operators, the solution would still 
have the useful property of error detection. In this case we can view the result in Eq. ( |28| ) 
as preservation for the part where the control bit is unity, and with error detection in the 
cases in which measurement of the control bit gives zero. When we measure zero for the 
control bit, we must discard the information in the channel. But even a single redundant 
parallel channel with the same input data would reduce the probability of information loss 
to the level of order (Tth) 2 - Since the latter is the essential level of error for the case in 
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which measurement of the control bit gives unity, this is the best result obtainable using our 
kind of mechanism. 

In our numerical example, we simply solve Eq. ( |25| ) for the density matrix evolution in 
the presence of noise, beginning at time t = and ending at time t = 2t e + t^, without 
breaking into the separate factors of the discussion above. This then includes every piece 
of Eq. (^), but with the additional corrections due to: a) Small non-adiabaticity (hopping) 
due to the fact that time taken by the encoding process is not infinite, b) Errors that occur 
from the noise interaction during encoding-decoding, due to the fact that the encoding is 
not instantaneous, c) The terms of higher order in Tt^, which matter when the probability 
of more than one error during the time th becomes appreciable. For our example, we choose 
the initial state for the information bit to be simply |1). Figure ^ shows the probability 
of a bit flip error arising from an error operator ( = a x ® 1 plotted as a function of the 
average number of applications of the error, Tt^, and on the same plot the results for the 
case in which the information is not encoded. We have also checked numerically to see that 
in cases in which we begin with a superposition of |1) and |0) for the information bit, and 
fine tune td to an appropriate value, the entire state, including the relative phase angle in 
the superposition, is preserved to the same degree as the bit-error values plotted in Fig. ||. 

Error prob. 

. 1 

.01 
. 001 
.0001 



0.01 0.1 1 * 

FIG. 6. The solid curve gives the bitffip error probability as a function of Tth for the case in 

which the encoding and decoding mechanisms are in place. The encoding time was taken to be 

t e = 0.01 x i/j. The dashed curve gives the error generated with the same noise source but without 

the encoding. 




IV. CONCLUSION 

We have considered two related topics in the study of closed and open quantum system 
evolution that go beyond one-body system Hamiltonians by including interaction terms be- 
tween the system particles. In the first topic, we have shown that entangled states that are 
symmetric under the permutation of the basis states can be generated via adiabatic tran- 
sitions from initially factorisable states which have this symmetry. In a two-qubit system, 
the outcome is the symmetrical Bell state -j= (| 10) + 1 01 ) ) , while the three-qubit analogue 
is the set of two W states. Noise in moderate quantitites generally reduces entanglement 
over time, and for such cases the limit of large noise sees entanglement remaining zero at all 
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times. However, for the special case of "side-blind" noise that leads to motional narrowing, 
the creation of entanglement is preserved despite a large system-environment interaction 
rate. 

In our second topic, we have shown that for small noise levels, we can use an entangling- 
pause-unentangling cycle to store and recover quantum data in a two qubit system, with 
error probability quadratic rather than linear in the parameter FT where F is the error rate 
for unencoded data and T is the time elapsed. Calculations of encoding that utilise the 
entanglement give a concrete demonstration of the compatibility of two requirements on the 
encoding time: a) that it be long enough achieve adiabaticity; b) that it be short enough to 
limit the effects of noise during the encoding process. 
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